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Correlations of two charged identical kaons (K ch K ch ) are measured in pp collisions at y/s = 7 TeV 
by the ALICE experiment at the Large Hadron Collider (LHC). One-dimensional K ch K ch correlation 
functions are constructed in three multiplicity and four transverse momentum ranges. The K ch K ch 
femtoscopic source parameters R and A are extracted. The K ch K ch correlations show a slight in- 
crease of femtoscopic radii with increasing multiplicity and a slight decrease of radii with increasing 
!><! transverse momentum. These trends are similar to the ones observed for kk and K^K^ 1 correlations 

^ ' m pp and heavy-ion collisions. However at high multiplicities, there is an indication that the one- 

dimensional correlation radii for charged kaons are larger than those for pions in contrast to what 
was observed in heavy-ion collisions at RHIC 
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1 Introduction. 

Extremely high energy densities achieved in heavy-ion collisions at the Large Hadron Collider (LHC) 
may entail the formation of the Quark-Gluon Plasma (QGP), a state characterized by partonic degrees 
of freedom flX Studying the QGP is the main goal of the ALICE experiment (A Large Ion Collider 
Experiment) 0. The system created in ultrarelativistic pp collisions at LHC energies might be simi- 
lar to the system created in non-central heavy-ion collisions because of the large energy deposited in 
the overlapping region and therefore may also manifest a collective behavior. The highly compressed 
strongly-interacting system is expected to undergo longitudinal and transverse expansions. Experimen- 
tally, the expansion and the spatial extent at decoupling are observable via Bose-Einstein correlations. 

Bose-Einstein correlations of two identical pions at low relative momenta were first shown to be sensitive 
to the spatial scale of the emitting source by G. Goldhaber, S. Goldhaber, W. Lee and A. Pais 50 years 
ago 0. The correlation method since developed and known at present as "correlation femtoscopy" 
was successfully applied to the measurement of the space-time characteristics of particle production 
processes in high energy collisions, especially in heavy-ion collisions (see, e.g. El[5j|6l)- Bose-Einstein 
correlations of identical particles were widely studied in heavy-ion collisions at the Relativistic Heavy- 
Ion Collider (RHIC) Q, and were found to confirm the hydrodynamic type of collective expansion of 
the fireball created in such collisions. In heavy-ion collisions the decrease of the correlation radii with 
increasing particle momentum was usually considered as a manifestation of a collective behavior of the 
matter created in such collisions J6). Event multiplicities reached in 7 TeV pp collisions at the LHC are 
comparable with those measured in peripheral A+A collisions at RHIC, making the study of the particle 
momentum dependence of the correlation radii an important test of the collectivity in pp collisions. 

The ALICE Collaboration has already studied two-pion correlation radii in pp collisions at 900 GeV 
|H and 7 TeV |9|. and K°K° correlation radii in pp collisions at 7 TeV 11101 . Two-pion Bose-Einstein 
correlations in pp collisions at yfs = 900 GeV and 7 TeV have been successfully described within the 
EPOS+hydro model ifTTTl . It was shown that the hydrodynamic expansion substantially modifies the 
source evolution compared to the "classical" EPOS scenario with independent decay of flux-tube strings, 
allowing one to describe the transverse momentum dependence of the correlation radii at high multiplic- 
ities. 

The main motivations for carrying out the present K ch K ch (K + K + + K~Kr) femtoscopy analysis are: 1) 
study the transverse mass, mj, dependence of the correlation radii ("wiT-scaling" is expected to be an 
additional confirmation of the hydrodynamic type of expansion J6l), 2) get a clearer signal (kaons are 
less affected by the decay of resonances than pions). 

Previous K ch K ch studies carried out in Pb-Pb collisions at SPS by the NA44 and NA49 Collaborations 
lfT2l and in Au-Au collisions at RHIC by the PHENIX Collaboration |[T3l revealed scaling in transverse 
mass: the source sizes versus mj for different particle types (n, K) fall on the same curve. 

K ch K ch studies were performed in a combined data from aa, pp, and pp collisions at ISR by AFS 
Collaboration [14], in e+e~ collisions at LEP by the OPAL and DELPHI Collaborations H21, in ep 
collisions by the ZEUS Collaboration lfT6l . Due to statistics limitations, only one-dimensional radii were 
extracted in these experiments, no multiplicity and transverse momentum studies were performed. 

In this article we present the measurements of Bose-Einstein correlations for charged kaons in pp colli- 
sions at y/s = 7 TeV performed by the ALICE Collaboration at the LHC. The present study is the first 
femtoscopic K ch K ch study to be carried out in pp collisions and in more than one multiplicity and pair 
transverse momentum, kj, range. 

The paper is organized as follows: in Section II we describe the ALICE experimental setup and data 
taking conditions for the data sample used in this work. In Section III we present the correlation mea- 
surements and the correlation functions. In Section IV we show the main results obtained in this work: 



2 



The ALICE Collaboration 



the one-dimensional radii extracted from the data. We discuss various observed features and compare the 
results with and %% radii previously measured by the ALICE Collaboration. Finally in Section V 
we summarize our results. 

2 Data analysis. 

Approximately 300 million minimum-bias events at yfs = 7 TeV, recorded in 2010, were analyzed. The 
ALICE Time Projection Chamber (TPC) and Inner Tracking System (ITS) were used for charged particle 
track reconstruction and the determination of the primary vertex of the collision. 

The TPC identifies charged particles according to their ionization trajectories in the Ne — CO2 gas. The 
ionization electrons drift up to 2.5 m from the central electrode to the end caps to be measured on 159 
padrows, grouped into 18 sectors. The position at which the track crosses the padrow is determined with 
a resolution of 2 mm and 3 mm in the drift and transverse directions, respectively. The ITS consists of 
six silicon layers, two innermost Silicon Pixel Detector (SPD) layers, two Silicon Drift Detector (SDD) 
layers, and two outer Silicon Strip Detector (SSD) layers, which provide up to six space points for each 
track. 

The forward scintillator detectors VZERO were included in the minimum-bias trigger and their timing 
signal was used to reject the beam-gas and beam-halo collisions. The minimum-bias trigger required a 
hit in one of the VZERO counters or in one of the two inner layers of the SPD. The VZERO detectors are 
placed along the beam line at +3 m and -0.9 m from the nominal interaction point. They cover a region 
2.8 < T] < 5.1 and -3.7 < r\ < -1.7. 

ALICE provides excellent particle identification capabilities, using the measurement of specific particle 
energy loss (dE/dx) in the TPC and the ITS and the time-of-flight (?tof) 

information obtained in the Time-Of-Flight (TOF) detector. The TOF detector is based on Multi-gap 
Resistive Plate Chambers (MRPCs) in a cylindrical configuration at radius 370-399 cm from the beam 
axis, with about 153.000 readout channels of dimension 3.5 x 2.5 cm 2 . The start time of the collision 
(event time zero) is measured by the TO detector, an array of Cherenkov counters located at +350 cm and 
-70 cm along the beam-line. If the TO signal is absent the start time is estimated from the particle arrival 
times at the TOF. The overall time-of-flight resolution depends on the TOF timing signal resolution 
(better then 100 ps), the accuracy of the reconstructed flight path and the uncertainty in the event start 
time. The resulting time-of-flight resolution is about 160 ps. 

The following criteria were applied for the event selection: 

- z-position of the reconstructed vertex within ±10 cm around the geometrical center of the ALICE 
detector; 

- at least one particle in the event reconstructed and identified as a charged kaon (In fact, the cor- 
relation signal is constructed from events having at least two same-charged kaons (a pair). The 
one-kaon events do contribute to the mixed background. It was verified that including one-kaon 
event to the mixed background does not change the shape of the correlation function. ) 

The criteria for track selection are listed below: 

- the kaons were selected in the kinematic ranges: |tj | < 1.0 and 0.15 < pj < 1.2 GeV/c; 

- tracks must include at least 70 space points (or clusters) out of a maximum possible number 159 
in the TPC, and two space points in the ITS (of maximum 6). 
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- the quality of a track was determined by the % 2 /N value for the Kalman fit to the constructed 
position of the TPC clusters (N is the number of clusters associated with the track); the track was 
rejected if the value was larger than 4.0 (2 degrees of freedom per cluster); 

- in order to reduce the number of secondary particles it was required that the particle trajectory 
distance from the primary vertex was less than 0.2 cm in the transverse plane and less than 0.25 cm 
in the beam direction. 

Usually the femtoscopic correlation functions of identical particles are very sensitive to two-track re- 
construction effects because particles have close momenta and thus close trajectories. The "splitting" of 
the tracks means that one track was reconstructed as two, and "merging" means that two different tracks 
were reconstructed as one. For the correlation structures measured in pp collisions, with characteristic 
widths ~ 0.2 GeV/c, track splitting and track merging in the event reconstruction are small effects, but 
we applied the standard femtoscopic double track cuts (see for details ll9l): 

- "anti-splitting cut": pairs which share more than 5% of clusters in the TPC were removed; 

- "anti-merging cut": pairs that are separated by less than 3 cm at the entrance of the TPC were 
removed. 

Pair cuts were applied in exactly the same way for real (signal) and mixed (background) pairs. 

In the present analysis the limit pj < 1.2 GeV/c for kaon selection on the TPC and TOF signals was 
used in order to ensure a high purity of the kaon sample. Kaons were selected by requiring that the 
deviation of the specific (dE/dx) energy loss in the TPC from that calculated with a parametrized Bethe- 
Bloch formula be within some number of sigma standard deviations (N ajpc ). A similar A^tof method was 
applied for the particle identification in the TOF using the difference between the measured time-of-flight 
and the calculated one as a function of the track length and the particle momentum at each tracking step 
and for each particle mass. 

More details on the particle identification are given in Ref. ifTTl . where it is shown in particular that the 
fraction efficiency of the particles reconstructed by the TPC with associated signal in the TOF (the TOF 
matching efficiency) deviates from 30% up to 55% in the /?T-region - under study in our paper. In the 
present analysis strict cuts on TPC and TOF signals for kaon selection were used in order to provide the 
better purity of the kaon sample. The relative contribution of the kaons from the sample with TOF signal 
to the full sample of the identified kaons used in this analysis is about 60%. 

In the present analysis strict cuts on TPC and TOF signals for kaon selection were used in order to 
provide the better purity of the kaon sample. If the TOF signal was not available, the following cuts were 
taken: 

- N ajpc < 1 for p < 0.35 GeV/c; 

- N ajpc < 2 for 0.35 < p < 0.6 GeV/c; 

- at larger momenta the tracks were rejected because of significant pion contamination. 
If the TOF signal was available, we required that N aj0P < 3 and < 3. 

Figure Q] shows the transverse momentum dependence of the kaon purity (the ratio of correctly identified 
kaons to all identified ones) obtained with Monte Carlo (PYTHIA) simulations. The contamination 
comes mainly from e + /e~ with maximum ~ 25% for 0.35 < pi < 0.6 GeV/c and also from pions at 
the level ~ 1 — 3% for 0.35 < pj < 1.2 GeV/c. In consequence, the probability of selecting an ee or 
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Fig. 1: Purity and contaminations of selected kaons in pp collisions at s/s = 7 TeV. 

a %% pair instead of a KK pair, even in the case where the contamination is maximal, is still rather 
low, smaller than ~ 6%. The probability of selecting an eK(Ke) pair instead of a KK pair can reach 
~ 37% for 0.35 < pi < 0.6 GeV/c. The probability of selecting a 7rK(K7r) pair instead of a KK pair 
is less then ~ 6% for 0.35 < p? < 1.2 GeV/c. Note that such contamination modifies only the strength 
of correlation and not the shape of the correlation function. Using the Monte Carlo information and 
also the purity of kaons estimated from the measured TPC dE/dx distributions and a realistic single- 
particle kinematics we have estimated the purity of kaon pairs at low relative momentum as: p = 
0.84±0.05 (stat.) ±0.15 (syst.), 0.61 ±0.03 (stat.) ±0.12 (syst.),0.79±0.04 (stat.)±0.07 (syst.),1.0± 
0.05 (stat.) ±0.05 (syst.) for pair transverse momentum k T = |p TJ +p Ti2 |/2: (0.2-0.35), (0.35-0.5), (0.5- 
0.7), (0.7-1.0) GeV/c, respectively. 

3 Charged kaon correlation functions. 

Momentum correlations are usually studied by means of correlation functions of two or more particles. 
Specifically, the two-particle correlation function CF (pi,p2) =A(pi,p2) /S(pi,P2) is defined as the ratio 
of the two-particle distribution in a given event A(pi,p 2 ) to the reference one, 5(pi,p 2 ), where pi and 
P2 are the momentum vectors of the two particles. In the present analysis the reference distribution is 
constructed by mixing particles of a given class, as described below. 

The analysis was performed in three multiplicity ranges based on the measured charged-particle multi- 
plicity, N c h- (1 — 11), (12 — 22), (> 22), and in 4 ranges of pair transverse momentum kj: (0.2-0.35), 
(0.35-0.5), (0.5-0.7), (0.7-1.0) GeV/c. Event multiplicity was determined as the number of charged parti- 
cles emitted into the pseudorapidty range \r\ | < 1 and transverse momentum range 0. 12 < pj < 10 GeV/c. 
For each class of events we calculated the charged-particle pseudorapidity density dA^h/dr] corrected for 
the detection efficiency obtained with Monte Carlo. The considered event multiplicity ranges (1 — 11), 
(12 — 22), > 22, correspond to mean charged particle densities, d/Y c h/dT], of 3.2, 8.1 and 17.2, respec- 
tively, with systematic uncertainties of ~ 5%. 

The numerators and denominators of positive and negative kaon distributions were summed up before 
constructing the ratio (K + K + and K~K~ correlation functions were found to coincide within errors thus 
justifying the procedure). The function is normalized to unity in the range 0.5 < q mw < 1.0 GeV/c, where 
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<7inv = y\q\ 2 -ql, q = pi- P2, and qo = E\ — E^. The range for normalization was chosen outside the 
Bose-Einstein peak. 

The correlation function is fitted by a single-Gaussian |[T8l : 

CF(q im ) = (1 - A + K{q im ) (A exp (-tfjUL) ) ) D{q inv ), (1) 

where the factor K(q mv ) is the Coulomb function integrated over a spherical source of 1 fm. The function 
D(qinv), "baseline", takes into account all non-femtoscopic correlations, including the long-range corre- 
lations due to energy-momentum conservation. The parameters R[ m and A describe the size of the kaon 
source, and the correlation strength, respectively. The A -parameter depends also on purity and decreases 
if the purity is not 100%. The R\ m is measured in the pair rest frame. 

The baseline was fitted by a standard quadratic polynomial 

£K<7inv) = 1 + aq™ + bqLv (2) 

To estimate the systematic errors due to the fitting procedure, other functions with derivatives equal to 
zero at q mv = were also employed, such as: 

D (*nv) = y/l+aq^ + b^ (3) 

and the Gaussian: 

D(q- mv ) = b(l+exp(-aql w )). (4) 

The PERUGIA-2011 tune jTJij of the Monte Carlo event generator PYTHIA (20l describes well the 
kaon spectra in pp collisions at LHC energies. Therefore, it was used to simulate the correlation function 
without the Bose-Einstein effect. 



4 Results and Discussion. 

Figure [2]presents the experimental two-kaon correlation functions and those obtained from a simulation 
using PERUGIA-201 1 (open circles) as a function of the invariant pair relative momentum. As one can 
see, the Monte Carlo simulation reproduces well the experimental correlation function at large <7i nv , i.e. 
the long-range correlations. The model does not contain the Bose-Einstein effect, so the enhancement 
at low q mv is due to non-femtoscopic correlations in PYTHIA, probably arising from mini-jets. The 
baseline points, obtained from PERUGIA-201 1 were fitted to Eq. ©. The parameters a and b were used 
in the fitting of the experimental points by Eq. £T|). The same method was used to model the baseline for 
the ALICE 7171 correlation studies in 0.9 [81 and 7 TeV [9| pp collisions. Figures [3]|4] and Table [Qpresent 

the one-dimensional A -parameters and Gaussian radii versus mj = \Jk^ + m\ including statistical and 
systematical errors. 

In order to estimate the systematic error from the choice of baseline functional form we repeated the 
fitting procedure using the baseline fitted with Eqs. (j3]|4). The radii obtained in the three ways differ 
by less than 4 % at low multiplicities (or kj) and by up to 10% at high multiplicities (or kj). The 
systematic errors estimated from varying the q mv fit range are below 2 % and up to 15% at low and high 
multiplicities (or kj) bins respectively. The systematic errors from splitting and merging effects were 
estimated by using different start points for the fit of correlation function: 0.03 and 0.06 GeV/c, and its 
are about 2-6%. The systematic error connected to the Coulomb function in Eq. ([TJ) was calculated in 
the following way: at first, the radius of the spherical source was taken equal to 1 fm, then the fitting 
procedure was repeated using these radii ±3c# (where Or is the total error) as the argument of the 
Coulomb function. The obtained systematic error is about 2-4%. 
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Fig. 2: Correlation functions versus qi nv for identical kaons from pp collisions at y/s = 7 TeV (solid circles) and 
those obtained with PERUGIA-2011 (open circles). Positive and negative kaon pairs are combined. The three 
columns represent the samples with different charged-particle multiplicities: (1 — 11), (12 — 22), (> 22), the four 
rows represent the four pair transverse momentum ranges: (0.2-0.35), (0.35-0.5), (0.5-0.7), (0.7-1.0) GeV/c. The 
lines going through the points represent the Gaussian fits discussed in the text. 
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Fig. 3: A-parameters of K ch K ch versus mj extracted by fitting correlation functions shown in Fig. |2]to Eq. (Q} and 
the baseline to Eq. ©. For comparison the K^K^ [10| A-parameters measured by ALICE in 7 TeV pp collisions 
are also shown. Statistical (darker lines) and total errors are shown. The points corresponding to the second and 
third multiplicity bins are offset by 0.03 GeV/c 2 for clarity. 



The mj dependence of A shown in Fig.[3]demonstrates that A varies within the range ~ 0.3 — 0.5 (except 
the first point at lowest multiplicity and lowest kj, which is ~ 0.2). As seen in Fig.[3j the A-parameters for 
K ch K ch are generally smaller than those for K°K° ifTOl . There are several reason for the A-parameter to be 
less than the ideal case of unity. Possible causes may be such as a partially coherent source, a contribution 
to the observed kaons from decays of long-lived resonances, a deviation from the Gaussian parameteri- 
zation due to a mixture of sources with different radii (see e.g. J2TJ), or a particle misidentification. The 
latter cause influences mostly the K ch K ch sample in the momentum range 0.4 < pj < 0.6 GeV/c (Fig.Q]). 
For the future, it seems desirable to improve the PID and thus increase the purity and to improve the 
description of non-femtoscopic correlations using different models. 

The K ch K ch correlation radii in Fig.@]show an increase with multiplicity in agreement with the KK radii at 
900 GeV El, and 7 TeV El, and the K°K° radii ED, as it was observed for nn correlations in heavy-ion 
collisions |f22| . These radii also decrease with increasing mj for the large multiplicity bins A c h (12 — 22) 
and A c h (> 22). Such a tendency was found for pions [9] and neutral kaons K ( S 'K° iflOl in pp collisions, 
and pions in heavy-ion collisions at LHC energies Il22l . In the low multiplicity bin A c h (1 — 11) charged 
kaons show a completely different ^-dependence of the radii: these radii increase with kj. This effect 
is qualitatively similar to that of pions 

It was observed that A-parameters (Fig. [3) are correlated with the radii (Fig. Such a correlation can 
result from the non-perfect fit results if the fit quality depends on mj. The reasons of such a dependence 
are: non ideal description of the baseline by PERUGIA-201 1, especially at large kj, non-Gaussian shape 
of the source due to resonance contribution and non-spherical shape of the source. The last two points 
mean that our one-dimensional Gaussian fit is only an approximate description of the source. In pp 
collisions the effect of this non-Gaussian shape of the correlation function due to different sizes in the 
'x-y-z' directions plays a more important role than in heavy-ion collisions. This requires a detailed 3D 
analysis, which is foreseen for K ch K ch correlation functions with the new large set of data recorded by 
the ALICE Collaboration in 2011 and 2012. 
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Fig. 4: One-dimensional charged kaon radii versus ntj extracted by fitting correlation functions shown in Fig. |2] 
to Eq. (Q]) and the baseline to Eq. (0. For comparison the and K^K^ ITOl radii measured by ALICE in 7 TeV 

pp collisions are also shown. Statistical (darker lines) and total errors are shown. The points corresponding to the 
second and third multiplicity bins are offset by 0.03 GeV/c 2 for clarity. 



The mj dependence of the radii in heavy-ion collisions was inteipreted as the manifestation of the strong 
collective hydrodynamic expansion of the created matter [6]. The observed similar behavior in pp col- 
lisions, shown in FigJH has some specific features: 1) at low multiplicity the radii increase with kj, 2) 
there is no distinct mj scaling: the kaon radii seem to be larger than the pion ones. The model calcula- 
tions performed in ifTTI can successfully describe the different behavior of pion correlation radii in low 
and high multiplicity bins, suggesting that the contribution of the hydrodynamic phase is negligible in 
low-multiplicity events, while for events with high multiplicity, it is substantial. 

As shown in l23l . due to the small size of the created system in pp collisions, the flow of resonances 
may play a significant role in large multiplicity bins, where essential hydrodynamic collective flow is 
expected ifTTIl . According to simple chemical model calculations ifTOl . the influence of this flow should 
be relatively smaller for kaons than for pions, leading to the effect that the kaon radii can be larger than 
the pion ones. The measured K ch K ch correlation radii displayed in Figj4] support such an hypothesis, 
however a detailed theoretical study is needed. 



5 Summary. 

The ALICE Collaboration has measured charged kaon correlation functions in pp collisions at yfs = 
7 TeV at the LHC. In agreement with the previous measurements in pp and heavy-ion collisions at lower 
energies, the extracted correlation radii R mv increase with the event multiplicity and decrease with the 
pair transverse mass/momentum. The novel features are some hints to the increase of the radii with mj in 
the low-multiplicity bin and to the fact that kaon radii are larger than the pion ones. These peculiarities 
deserve further experimental and theoretical studies. 



Charged kaon femtoscopic correlations in pp collisions at = 7 TeV 



9 



Table 1: K ch K ch source parameters vs. kj for \fs—l TeV pp collisions. Statistical and systematic errors are listed. 
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